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Abstract 

We study the radial Schroedinger equation for a particle of mass m in the field of a singular attrac- 
tive a/r^ potential with 2ma > 1/4. This potential is relevant to the fabrication of nanoscale atom 
optical devices, is said to be the potential describing the dipole-bound anions of polar molecules, 
and is the effective potential underlying the universal behavior of three-body systems in nuclear 
physics and atomic physics, including aspects of Bose-Einstein condensates, first described by Efi- 
mov. New results in three-body physical systems motivate the present investigation. Using the 
regularization method of Beane et al., we show that the corresponding "renormalization group 
flow" equation can be solved analytically. We find that it exhibits a limit cycle behavior and 
has infinitely many branches. We show that a physical meaning for self-adjoint extensions of the 
Hamiltonian arises naturally in this framework. 

PACS numbers: 34.20.-b, 03.65. Ca, 03.65.Ge, ll.lO.Hi, 
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In this note, we study the regularization and renormahzation of the singular attrac- 
tive ajr^ potential, a problem motivated, in part, by recent results in three-body physics. 
Specifically, there has been renewed interest in the non-relativistic three-body system with 
short-range interactions. The current investigations are stimulated by the promise that ef- 
fective field theories (EFT'S) offer a systematic and model-independent treatment of atomic, 
nuclear and hadronic physics at low energies. That is, a low-energy system with a clear- 
cut separation of distance scales can be described by an EFT involving explicitly only the 
long-wavelength degrees of freedom. The short-range dynamics can be treated as a set of 
local operators which correspond to delta-function interactions in coordinate space. The 
details of the short distance physics cannot be of importance to the low energy aspects of 
the system; if they are, it is an indication of a need for renormahzation of the EFT. 

Such a renormahzation of an EFT of a three-body system with delta-function two-body 
interactions |])] has lead to the rediscovery of the one-parameter contact three-body interac- 
tion shown to restore the lower bound of this three-body Hamiltonian j^. (Here we mention 
that the unboundedness of the Hamiltonian from below is interpreted in EFT's as the onset 
of short-distance physics whose effect must be included in local counterterms jsl; hence the 
three-body counterterm in the EFT equations for 3-body systems, both bound and scatter- 
ing states [41.) The three-body counterterm exhibits a convergence of the renormahzation 
group fiow to one-dimensional limit cycles. This was earlier proved in Ref. Q], a mathemat- 
ical analysis of the Efimov effect which occurs in bound 3-body systems when more than 
one of the two-body subsystems has a zero energy resonance 1^1. The presumed novelty of 
a renormahzation group fiow with a limit cycle has inspired additional recent work 0]. 

As the applications of EFT's continue 0, Q|, it will be important to understand how 
to explicitly renormalize higher orders in an EFT. For example, a renormalized equation 
for two-nucleon systems with explicit pion-exchange would be of great potential value. Pion 
exchange gives rise to a singular potential and the questions arise: Can the resummation 
of pion graphs be renormalized by a single local operator? Would this operator exhibit a limit 
cycle as does the three-body contact operator in the pion-less three-nucleon EFT? Could 
one calculate the evolution of such an operator analytically? A positive answer to the last 
question would help future numerical work with EFT's. Already, the short distance physics 
of the ^5*1 coupled channels of the single pion exchange potential has been renormalized 
by a short range four-nucleon counterterm using the method introduced in ref. [3], but the 
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treatment was numerical and these questions were not addressed in that investigation 

Such questions, coupled with the interesting limit cycle behavior found in the three-body 
system, have prompted investigations of the renormalization group behavior of the short 
range counterterms which serve to regularize given long range potentials (including singular 
potentials) in the two-body Schroedinger equation. Since these long range potentials are 
often singular at the origin, it has been argued that the short range interaction should not 
be represented by a 3-dimensional delta-function at the origin. Birse et al. |9] choose a delta- 
shell potential and Beane et al. [sj suggest that a simple attractive square well represents 
a "smeared out" delta-function potential. In either case, it is said, the details of the short 
distance regularization should not matter; the low energy aspects of the system should be 
invariant in the same way under suitable changes of the short range potential. 

For completeness, we list other regularization and renormalization schemes which do not 
follow from the separation of scales of an ETF, but also have been applied to the inverse 
square potential description of physical systems. The problem of a neutral atom interacting 
with a charged wire |ly| , relevant to the development of nanoscale atom optical devices, has 
been treated with the method of self-adjoint extensions 1J|. A short distance cutoff scheme 
which renormalizes the strength of the potential, yields a critical dipole moment that 
has been confronted with the experimental capture of electrons by dipole molecules and 
formation of anions as an example of quantum mechanical symmetry breaking 
These alternative regularizations are not the subject of our present investigation. 

In this note, we follow the regularization method of Beane et al. to obtain analytically 
the renormalization group behavior of the coupling constant of the short range attractive 
square well they use to regularize the long range inverse square potential. The potential 
is on the boundary between singular and regular potentials and thus does or does not require 
a self-adjoint extension, depending on the strength of the interaction. More interesting from 
the EFT point of view is the Efimov observation ^ that the low energy behavior of three- 
body systems is determined by a long range three-body effective interaction of the form 
1/ where R is built from the relative distances between the particles. Thus, the two-body 
inverse square potential is the analogue to the interaction in a three-body system in the 
limit of zero energy resonance (and infinite two-body scattering lengths). 

A further advantage of our analytic approach to the "EFT style" renormalization of the 
inverse square potential is that we can then more readily make contact with the mathe- 
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matically rigorous and well studied approach to regularization via self-adjoint extensions. 



The theory of self-adjoint extensions o: 
cycle behavior in three-body systems 



Hamiltonians underlies the first discussions of limit 



The self-adjoint extensions of the inverse square 



potential are well known [15[ and can be compared with the results of our study, which we 
now begin. 

The starting point of our study is the s— wave reduced radial Schroedinger equation for 
one particle of mass m in the external potential V{r): 



— - 2mV{r) + /c^ U = (1) 



where V{r) is given by 



3: 



dr 



V{r) = - - ' - > K,a>0). (2) 

That is, the long range attractive a/r^ second term in eq. (2) is cutoff at a short distance 



radius R by an attractive square well. As in we first solve eq. (1) for the zero energy 

solution {k = 0) ipo- It is given by: 

a/2 



ipo{r) = A— jj^ cos In— + (poj r>R (3) 

cos(i/ln-^ + /i?i/2\ 

where u = {2ma — 1/4)"^''^, 0o is the zero energy phase = {2mas)/R^ and Vo is an 

arbitrary scale. 

The usual matching condition of the wave function and its derivative at r = i? then 
yields: 

(2mas)"^^^ cot{(2mas)"^''^} = z/tan (uln (—^ + (po) (5) 

Following we now consider eq. (0) to be a transcendental equation defining the value of 
the short range coupling constant a^. This equation is of the form f3cot(3 = l/u and can 
be solved exactly in closed form usin g a method based upon the solution to the Riemann 
problem in complex variable analysis Iq . 

The solution to eq.© then turns out to be jl?! : 

Po = ±- — exp - / arg A, t - ] , u>0 6 

u \n Jo t 



(l dt\ 
invrexp — / axgVtnif) — , — oo < c^; < -|-oo, n = 1, 2, ... (7) 
yvr JO t J 



where: 



Ao(t) 



1 




{2mas) 



(8) 
(9) 
(10) 



\{t) 



1 + 



(11) 
(12) 



Formula (9) does not restrict uj to be positive, so that all the solutions that we consider follow 
from the positive branch of (7). If we interpret eq. (jSJ to be an equation for the running 



values of R/to except for those points where uj (or its inverse) vanish. This is potentially 
important when comparing our analytic solution with numerical solutions to eq. (5). Indeed, 
numerical investigations in general will mix solutions that only exist in a limited range of 
-R/ro-values {[3o in eq.(6)) with solutions (/?„ in eq.(7)) which we consider to be the only 
physically relevant solutions. From now on, we shall therefore use (3 to mean only a solution 
of eq. (7), and suppress the subscript n when it is not needed for the discussion. We then 
find from eq. (7) that eq. © has infinitely many roots, in agreement with 0]. We have 
plotted (3 in eq. © as a function of Inx (x = R/to) for fixed (po = 1-0, and = 1 in figures 
1 and 2. The strength of the long range potential a/r^ increases with succeeding figures; 
u = 0.5 in fig. 1, u = 3.0 in fig. 2. The running coupling j3 exhibits a limit cycle behavior 
for all values of u > 0; the period becomes smaller as the strength of the attractive a/r^ 
potential increases, according to the argument (z/lnx + 0o) of the tangent function in the 
source term of eq. (5). The behavior described by (3, for large enough u, is of a "sawtooth" 
type, with a periodic sharp increase of the value of the coupling constant with decreasing 
values of Inx. This "increase" is actually a genuine discontinuity of /3 at the zeros of 1/u; 
and must be a multiple of vr. Indeed, a discontinuity can only occur at a zero of cot (3 in 
order that (3 cot (3 be continuous at all points where l/uj is a continuous function. Once f3 
has reached the smallest positive zero {(3 = 7r/2), for some x— value, it increases by tt as x 
is further decreased. Altering the zero energy phase shift 0o from through 2 (and keeping 
n=l and z/ = 3) does not qualitatively change the appearance of the pattern of Fig. 2: the 
discontinuity moves to lower x , but the magnitude of the discontinuity remains the same 
multiple of n. This feature can be traced to the periodicity of the right hand side of eq. (0) 



coupling constant (3 = {2mas) 



'2 

, this amounts to requiring that f3 must be defined for all 
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with respect to 0o- The branches which correspond to roots with n > 1 are quahtatively 
similar to Fig 2., but the slow fall-off with decreasing x seen in Fig. 2 increases and the 
"saw-tooth" appearance becomes more of a rounded off square wave. This feature of the 
solutions is illustrated by Fig. 3 which plots (3 for the same strength and initial phase z/ = 3, 
0o = 1 as Fig. 2, but n has increased to 16. The magnitude of the discontinuity of this Piq 
remains vr, however. 

Finally we plot as Fig. 4 the analytical solution of eq. (0) for the values u = 2.0, 
0o = 0.0 to compare with the numerical solutions of eq. (0), with the same input, displayed 
in Figure 1 of Ref. j^. The latter numerical solution is in excellent agreement with the 
analytical solution presented here, if one allows a solution to go from a higher branch (n = 
2) to the next lower branch (n = 1) as it crosses a zero of 1/uj. However, as discussed 
above and shown on Fig.4, j3i itself must ultimately increase by tt after reaching its lowest 
value 71/2, thus exhibiting a limit cycle behavior of the solution.lt is evident that the limit 
cycle behavior of the solutions of equation (5) is a consequence of the requirement that a 
solution be defined for all values of R/vo for which omega or its inverse is nonzero. This 
requirement includes solutions for R ^ and therefore corresponds to our understanding 
of the emulation proposed in ref [3] of the contact term which encapsulates the short range 
dynamics of an EFT. From Figures 1-3, it is clear that the limit cycle behavior of a given 
branch continues to the left as i? ^ and Inx becomes arbitrarily small. Consider, however, 
the behavior for small i? of a numerical solution of ref. [3], shown in Figure 4 which segues 
smoothly between different branches of the analytic solutions of equation (5) as it crosses a 
zero of l/uj. For some value of negative In x, as R becomes arbitrarily small, the numerical 
solution must pass to the f^Q solution. But, as we have noted, this solution exists for only 
a limited range of R/vo values. Thus, the two numerical solutions of Fig. 4 would appear 
to not be defined for arbitrarily small R/vo- One could, however, choose another numerical 
solution corresponding to a higher value of n which does not pass to /?o on the way from 
large to arbitrarily small R/vo and avoid this problem. This exercise need not be performed, 
however, with our requirement of a well defined solution for all values of R/tq for which 
omega or its inverse is nonzero. 

A motivation fo. both the study of Ref. 3 and the present diseuss.on is the expeotat.on 
that the two-body inverse square potential is the analogue to the interaction in a three-body 
system in the limit of zero energy resonances. Indeed, the approximate solution of eq. (5) 
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FIG. 1: The running coupling constant /3 as a function of Inx = In ^ for (j)o = 1.0, u = 0.5, n = 1 
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FIG. 2: The running couphng constant /3 as a function of Inx for (po = 1-0, v = 3.0, n = 1 

displayed in eq. 8 of 3 is quite similar to the equation which describes the running of the 
three-body counterterm of the pion-less three-nucleon EFT's of Refs. 0,1^. Both equations 
have poles and the three-body counterterm seems to reach arbitrarily high values. We, to 
the contrary, find no poles in the analytic solutions of eq. (5). Furthermore, no evidence 
of multiple branches was found in the renormalized pion-less three-nucleon problem. The 
renormalization of short distance physics in these two problems needs more understanding 
in light of the results of Efimov . 

Now we turn to the bound state aspects of the related two problems (three-body system 
with contact potentials and the singular potential) and again find discrepancies. First 
we show that the regularization method that was used to solve the Schroedinger equation 
(1) with potential (2) amounts to specifying a particular self-adjoint extension in Case's 



solution of the bound state (B.S.) spectrum of the attractive singular potential [l5|, isl ]. 
In order to do this, we note that the B.S. wavefunction that solves eqs. (l)-(2) is given by: 
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FIG. 3: /? as a function of Inx for <^o = 1.0, = 3.0, n = 16 
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FIG. 4: /3 as a function of Inx for 0o = 0.0, = 2.0. The three branches, from bottom to top, 
correspond to n=l (dash-dotted curve), n=2 (dotted curve) and n=3 (solid curve). As explained 
in the text, the discontinuity in /? is always tt for each curve. The two thicker curves are two 
numerical solutions for the same parameters taken from ref. [3], as discussed in the text. 



where: 



= Cr^/^Ki^{kr) r> R (13) 
ijj = C'sm(Kr) r<R (14) 



K^^^-^-k^ (15) 



and C and C are constants. 

1 /2 

For kR << 1, the matching condition now gives (one still has KR = {2'mas) in that 
limit) : 

,= -Lexp^» + "^r(l + |.)-(»+l/2). (16) 
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where we have used eq. (0) together with the small-r behavior of ilj{r) 



■?/'(r) ~ r^''^ sin ^z/log ^ argr(l + zz/)^ . (17) 

The spectrum given in eq. (|T6|l is essentially the spectrum given by Case Q], where 
Case's arbitrary phase (which fixes the self-adjoint extension) B is now given by: 

B = (f)o + aTgT{l + iiy) (18) 

It is important to note that the binding energy Eb = {k'^)/2m, after this regularization, no 
longer depends on the cut-off radius R (for kR << 1) but instead on the arbitrary scale Tq. 
Thus, fixing the zero energy phase of the wavefunction (po removes the cut-off dependence of 
the B.S. spectrum for kR << 1 . Our result explicitly shows that the physical interpretation 
of the phase characterizing the self-adjoint extension of the Hamiltonian indeed can be found 
from a renormalization of the short-range coupling constant in the regularization method 
described in Ref. jsj]. By the same token, it illustrates why a cut-off method with a constant 
strength fails to provide a physical meaning for this arbitrary phase. 

Note, however, that the ground state of the potential remains at negative infinity, 
and no renormalization of the bound state spectrum has been achieved. Contrast this result 

with the restoration of the lower bound of 



A clue to this discrepancy may 



of the EFT style renormalization of Ref. 

the pion-less three-body problem obtained in Refs. 

ie in the distinction between the contact interaction used in Ref. 2 and the "EFT" type 
1 

3|] regularization of the conventional 3-dimensional delta function. That is, the attractive 
square well in eq. (2) does not provide a unique way of regularizing ("smearing out") 
a 3-dimensional delta function, and its limit when i? ^ is not the contact interaction 
discussed in Ref. 2 and 14. In that respect, it would be quite interesting to reexamine the 
solution of the Schroedinger equation with a singular a/r"^ potential in conjunction with 
local realizations of the contact interaction of Ref. 14 implemented by Kruppa, Varga and 
Revai [2^. Such a study might throw additional light on the corresponding renormalization 
group fiow properties in the 3-body problem. 
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